home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dpoequ.z / dpoequ
Encoding:
Text File  |  2002-10-03  |  3.8 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPOOOOEEEEQQQQUUUU((((3333SSSS))))                                                          DDDDPPPPOOOOEEEEQQQQUUUU((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPOEQU - compute row and column scalings intended to equilibrate a
  10.      symmetric positive definite matrix A and reduce its condition number
  11.      (with respect to the two-norm)
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DPOEQU( N, A, LDA, S, SCOND, AMAX, INFO )
  15.  
  16.          INTEGER        INFO, LDA, N
  17.  
  18.          DOUBLE         PRECISION AMAX, SCOND
  19.  
  20.          DOUBLE         PRECISION A( LDA, * ), S( * )
  21.  
  22. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  23.      These routines are part of the SCSL Scientific Library and can be loaded
  24.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  25.      directs the linker to use the multi-processor version of the library.
  26.  
  27.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  28.      4 bytes (32 bits). Another version of SCSL is available in which integers
  29.      are 8 bytes (64 bits).  This version allows the user access to larger
  30.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  31.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  32.      only one of the two versions; 4-byte integer and 8-byte integer library
  33.      calls cannot be mixed.
  34.  
  35. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  36.      DPOEQU computes row and column scalings intended to equilibrate a
  37.      symmetric positive definite matrix A and reduce its condition number
  38.      (with respect to the two-norm). S contains the scale factors, S(i) =
  39.      1/sqrt(A(i,i)), chosen so that the scaled matrix B with elements B(i,j) =
  40.      S(i)*A(i,j)*S(j) has ones on the diagonal.  This choice of S puts the
  41.      condition number of B within a factor N of the smallest possible
  42.      condition number over all possible diagonal scalings.
  43.  
  44.  
  45. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  46.      N       (input) INTEGER
  47.              The order of the matrix A.  N >= 0.
  48.  
  49.      A       (input) DOUBLE PRECISION array, dimension (LDA,N)
  50.              The N-by-N symmetric positive definite matrix whose scaling
  51.              factors are to be computed.  Only the diagonal elements of A are
  52.              referenced.
  53.  
  54.      LDA     (input) INTEGER
  55.              The leading dimension of the array A.  LDA >= max(1,N).
  56.  
  57.      S       (output) DOUBLE PRECISION array, dimension (N)
  58.              If INFO = 0, S contains the scale factors for A.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDPPPPOOOOEEEEQQQQUUUU((((3333SSSS))))                                                          DDDDPPPPOOOOEEEEQQQQUUUU((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      SCOND   (output) DOUBLE PRECISION
  75.              If INFO = 0, S contains the ratio of the smallest S(i) to the
  76.              largest S(i).  If SCOND >= 0.1 and AMAX is neither too large nor
  77.              too small, it is not worth scaling by S.
  78.  
  79.      AMAX    (output) DOUBLE PRECISION
  80.              Absolute value of largest matrix element.  If AMAX is very close
  81.              to overflow or very close to underflow, the matrix should be
  82.              scaled.
  83.  
  84.      INFO    (output) INTEGER
  85.              = 0:  successful exit
  86.              < 0:  if INFO = -i, the i-th argument had an illegal value
  87.              > 0:  if INFO = i, the i-th diagonal element is nonpositive.
  88.  
  89. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  90.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  91.  
  92.      This man page is available only online.
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.